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Abstract 

With the help of the general theory of the Heun equation, this paper completes previous work 
by the authors and other groups on the explicit representation of the massive gravitino propagator 
in four-dimensional de Sitter space. As a result of our original contribution, all weight functions 
which multiply the geometric invariants in the gravitino propagator are expressed through Heun 
functions, and the resulting plots are displayed and discussed after resorting to a suitable truncation 
in the series expansion of the Heun function. It turns out that there exist two ranges of values 
of the independent variable in which the weight functions can be divided into dominating and 
sub-dominating family. 
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I. INTRODUCTION 



The investigation of Green functions has always been at the heart of important develop- 
ments in quantum field theory and quantum gravity On the other hand, in recent years, 
developments in cosmology and string theory have led to renewed interest in supergravity 
theories in anti-de Sitter 2] and de Sitter space 

Thus, in our recent paper [4|, we performed a two-component spinor analysis of geomet- 
ric invariants leading to the gravitino propagator in four- dimensional de Sitter spacetime, 
following the two-spinor language [5| pioneered by Penrose. In that paper we also wrote 
down all the 10 different weight functions multiplying the invariants which occur in the 
massive gravitino propagator, relying upon the work by Anguelova and Langfelder jc|. It 
was also found there, that algebraically one can write down 8 weight functions denoted by 
a, ft, 7, 5, e, 9, r, u> in terms of a pair denoted by (tt, k), in case of de Sitter space. Going one 
step further, we also expressed k in terms of tt and tt' where tt was defined in this fashion: 
ix(z) = \[ztt{z) and tt(z) satisfies the Heun differential equation 3, S], whose solutions, 
denoted by Heun(a, q; b, c, d, e; z), with properly defined arguments, have in general four 
singular points, i.e. z = 0, 1, a, oo. 

In this paper we have explicitly written down, first, the algebraic expression of all the 
9 weight functions k, a, ft, 7, 5, e, 9, r, u in terms of tt(z) and tt'(z), where z is defined as 
z = cos 2 7^, with fi(x,x') being the geodesic distance between x and x' as defined in J4]. 
Finally, we will draw a few two-dimensional plots of these weight functions and classify their 
parameter space with respect to z in the region of our choice. 

The plan of this paper is as follows. In Sec. II we set up all symbo.s, basicahy recalling 
all relevant definitions of use in this paper from our previous one [4j|. Sec. Ill contains the 
explicit massive spin-3 / 2 propagator in four dimensions with all the ten invariant structures 
properly defined, along with the multiplicative weight functions written in terms of the (tt, k) 
pair. In Sec. IV we give a crash course on Heun's differential equations and write down 
several properties of the Heun function before showing that fr(z) satisfies a Heun equation 
with properly defined arguments that we will list there. Sec. V and the appendix are 
devoted to build a dictionary of all the 9 weight functions k, a, ft, 7, 5, e, 9, r, uj written in 
terms of tx(z) and tt'(z) only, where prime denotes derivative with respect to z instead of 
being the derivative with respect to /x, the geodesic distance function. Then in the last Sec. 
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VI we display several two-dimensional plots showing the functional behavior of each of the 
10 weight functions with respect to z. These show that there exist two ranges of values of 
z in which the weight functions can be divided into dominating and sub-dominating family. 
Moreover, it appears helpful to have the result of a lengthy calculation completely worked 
out. Eventually, we give further details on the plots in the section devoted to concluding 
remarks. 

In light of recent mathematical developments in 9|, it might be possible to expand the 
Heun functions in the gravitino propagator as a combination of finitely or infinitely many 
hypergeometric functions, which in turn occur in the more familiar formulae for bosonic 



propagators in de Sitter space [lOj. Thus, our work might help relating fermionic and 
bosonic propagators in four- dimensional de Sitter space through special-function techniques, 
double-checking the expectations from supersymmetry. 

II. A REVIEW OF A FEW USEFUL DEFINITIONS 

It has been more than two decades since Allen and co-authors used intrinsic geometric 
objects to calculate correlation functions in maximally symmetric spaces; their results, here 



exploited, were presented in two papers ll|, |l2j. In this section we would like to review, 



first, the elementary maximally symmetric bi-tensors which have been discussed previously 
by Allen and Jacobson [111 ] . More recently, the calculation of the spinor parallel propagator 



has been carried out in arbitrary dimension 



3- 



A maximally symmetric space is a topological manifold of dimension n, with a metric 
which has the maximum number of global Killing vector fields. This type of space looks 
exactly the same in every direction and at every point. The simplest examples are flat space 
and sphere, each of which has |n(n + 1) independent Killing fields. For S n these generate 
all rotations, and for MJ 1 they include both rotations and translations. 

We consider a maximally symmetric space of dimension n with constant scalar curvature 
n(n — 1)/R 2 . For the space S n , the radius R is real and positive, whereas for the hyperbolic 
space H n , R = il with I positive, and in the flat case, MJ 1 , R = oo. If we further consider 
two points x and x', which can be connected uniquely by a geodesic, with fi(x, x') being 
the geodesic distance between x and x\ then n a (x, x') and n a ' (x, x') are the tangents to the 



3 



geodesic at x and x', and are given in terms of the geodesic distance as follows: 



n a {x, x') = V a fi(x,x') and n a t(x, x') = V a >fi(x, x'). (2.1) 

Furthermore, on denoting by g a b ,(x,x') the vector parallel propagator along the geodesic, 
one can then write n b ' = —g b ' a n a . Tensors that depend on two points, x and x\ are bitensors 
[l^ . They may carry unprimed or primed indices that live on the tangent space at x or x' . 
These geometric objects n a , n a ' and g a b , satisfy the following properties 



V a n b = A(g ab - n a n b ), (2.2a) 
\7 a n w = C(g ab i + n a n v ), (2.2b) 
^a9bd = -{A + C){g ah n c , + g ac m b ), (2.2c) 



where A and C are functions of the geodesic distance /i and are given by [TiJ 



A=-cot^ and C = l — -, (2.3) 

R R Rsm(fj,/RY K J 

for de Sitter spacetime and thus they satisfy the relations 

dA/dfi = -C 2 , dC/dn=-AC and C 2 - A 2 = 1/R 2 . (2.4) 
Last, with our convention the covariant gamma matrices satisfy the property 

{r",r"} = 2Ig» v . (2.5) 

n 

In our previous work [4J , we followed the conventions for two-component spinors, as well as 



all signature and curvature conventions, of Allen and Lutken [121 ]. and hence we used dotted 
and undotted spinors instead of the primed and unprimed ones of Penrose and Rindler |j] . 
In our work a primed index indicates instead that it lives in the tangent space at x', while 
the unprimed ones live at x. The fundamental object to deal with is the bispinor Dj±'(x, x') 
which parallel transports a two-component spinor (p A at the point x, along the geodesic to 
the point x' , yielding a new spinor x A ' a t x', i.e. 

X A ' = <P A Df(x,x'). (2.6) 

Complex conjugate spinors are similarly transported by the complex conjugate of Df (x, x'), 
which is (x,x'). A few elementary properties of Dj± were listed in Sec. IV of [4j. It 



is worth mentioning that the covariant derivatives of the spinor parallel propagator were 
defined to be 



V AA Di ' = (A + C) 



In . n b' _ . Tj B' 

2 AA B n BA LJ A 



(2.7) 



where A and C are defined in (2.3). 

The two basic massive two-point functions for spin- 1/2 particle, were defined by 



A A'B' 



P = {<j> A {x)f (x')} = f^)D A A ,n 



and in de Sitter space they turned out to be 4]: 



pAB> 



A A'B' 



lim f DS (Z + ie)D A ,n 



Qf*= lim g DS (Z + ie)D 



-r^AB' 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



where (F) stands for the Feynman Green functions with Jds and gos defined in this fashion: 

Ids = N DS (1 - Z)V 2 F(a, b; c; Z), (2.12) 

g DS = -iN DS 2- 3/2 m\R\Z l/2 F(a, b;c + l;Z). (2.13) 

Moreover, after doing some algebra one can rewrite the final answer for the constant N^s 
as 



N. 



DS 



—i 


Rm 


(1 — m 2 


R 2 ) 


8V27T 


R 


3 sinh 7r 


Rm\ 



(2.14) 



We also note that F(a, b;c; Z) and F(a,b;c + 1;Z) are two independent solutions of the 
Hypergeometric equation [15J, |l6j : 



H(a,b,c; Z)w(Z) = 0, 
where H(a, b, c) is the hypergeometric operator 

H(a, b, c; Z) = Z(l - Z)-^ + [c - (a + b + l)Z}-^ - ab. 



(2.15) 



(2.16) 



5 



III. MASSIVE SPIN-3/2 PROPAGATOR 



In this section we consider the propagator of the massive spin-3/2 field. Let us denote 
the gravitino field by In a maximally symmetric state | s) the propagator is 

SfJ(x,x') = (smx)^,(x')\s}. (3.1) 

The field equations imply that S satisfies 

{T»» x D p - mY^) V*/V = * ( ^V ,/ *V (3.2) 

A. The ten gravitino invariants 

It is very convenient to decompose the gravitino propagator in terms of independent 
structures constructed out of n M ,?v , g^i and A^, [4J. Thus, the propagator can be written 
in geometric way following Anguelova et al. [6| (see also 

S\v> a f3> = a(» £ Ai /AV + P(ji) n x n u ,A a p, + 7 (/i) g Xv >(n a Y° A) Q > 

+5(/i) n A 7v(n CT r CT A)V + e(^) n A (ryA)V + 0(» rv(r A A)> 
+r(/i) n A (n .r°T^A) a /3 / + rv(n ff PT A A)> 

+tt(^) (r A rvA)> + k(^) Kr<r A ivA)v (3.3) 

B. The weight functions multiplying the invariants 

A rather tedious but straightforward calculation gives a system of 10 equations for the 
10 coefficient functions a, k in (I3.3P as found in (see equations (3.6)-(3.15) in [6j). It was 
also found there that one can easily express the algebraic solutions for a, (3, 7, 5, e, 8, r, uj in 
terms of the (tt, k) pair in case of de Sitter space, i.e. (hereafter we set n = 4 in the general 



formulae of 



6] , since only in the four-dimensional case the two-component-spinor formalism 
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can be applied) 



CO 



T 

e 
a 



2mCK+ ((A + C) 2 -m 2 )7r 

(m 2 + R- 2 ) ' 
((A-C) 2 -m 2 )K-2mCn 

(m 2 + R- 2 ) ' 
2mCn + ((A + C) 2 — m 2 )7r 

(m 2 + R- 2 ) ' 
-{[{A - Cf + 2/R 2 ] + m 2 )K + 2mCir 
(m 2 + R~ 2 ) : 



— r — Air, 
(3 = 2co, 
7 = e — 2k, 
5 = 2e + i(K-6), 

where we have used the relation C 2 — A 2 = 1/R 2 . 

Furthermore, from (13. 4p we can immediately see that 



T = UJ 



and 



e + 6 



-2k. 



(3.4) 



(3.5) 



IV. HEUN'S DIFFERENTIAL EQUATION: A PRIMER 



The canonical form of the general Heun differential equation is given by ( 8|, 181 ]) 

72„. /. g 



^+(2 + + 

dz 2 \z z — 1 z — a J dz z\z — l)[z — a) 



dy a(3z-q 

+ — — -y = 



(4.1) 



The four regular singular points of the equation are located at z — 0, 1, a, oo. Here a G C, 
the location of the fourth singular point, is a parameter (a ^ 0, 1), and a, (3, 7, S, e 6 C are 
exponent-related parameters. 

The solution space of the Heun differential equation is specified uniquely by the following 
Riemann P-symbol: 

1 d 00 
P < a ; z } ■ (4.2) 
1-7 1-8 1-e (5 

This does not uniquely specify the equation and its solutions, since it omits the accessory 
parameter q e C. The exponents are constrained by 



a + /3~7-5-e + l = 0. 



(4.3) 



This is a special case of Fuchs's relation, according to which the sum of the 2n characteristic 
exponents of any second-order Fuchsian equation on CP 1 with n singular points must equal 
n — 2 (ISO. 



There are 2x4 = 8 local solutions of (4.1) in all: two per singular point. If 7 is not a 
nonpositive integer, the solution at z = belonging to the exponent zero will be analytic. 
When normalized to unity at z = 0, it is called the local Heun function, and is denoted 
Hl(a, q; a, 3, 7, 8; z) [8J. It is the sum of a Heun series, which converges in a neighborhood 



of z = 
integer. 



201 ] . In general, Hl(a, q; a, (3, 7, 8; t) is not defined when 7 is a nonpositive 



If e = and q = add, the Heun equation loses a singular point and becomes a hyperge- 
ometric equation. Similar losses occur if 8 — 0, q — ad, or 7 = 0, q = 0. This paper will 
exclude the case when the Heun equation has fewer than four singular points. The case, in 
which the solution of (4.1) can be reduced to quadratures, will also be ruled out. If ad = 
and q = 0, the Heun equation (4.1) is said to be trivial. Triviality implies that one of the 
exponents at z = 00 is zero (i.e., ad = 0), and is implied by absence of the singular point 
at z = 00 (i.e., ad = 0, a + (3 = 1, q = 0). 



A. Reducing Heun to hypergeometric 

The transformation to Heun (S)) or hypergeometric (h) of a linear second-order Fuchsian 
differential equation with singular points at z = 0, l,d, 00 (resp. z = 0,1, 00), and with 
arbitrary exponents, is accomplished by certain linear changes of the dependent variable, 
called F-homotopies (see [3] and 0, §A2 and Addendum, §1.8].) If an equation with 
singular points at z — 0,1, a, 00 has dependent variable u, carrying out the substitution 
u(z) = z~ p (z — l)~ a (z — a)~ T u(t) will convert the equation to a new one, with the exponents 
at z = 0,1, d reduced by p, a, r respectively, and those at z = 00 increased by p + a + r. 
By this technique, one exponent at each finite singular point can be shifted to zero. 

In fact, the Heun equation has a group of F-homotopic automorphisms isomorphic 
to (Z 2 ) 3 , since at each of z = 0, 1, a, the exponents 0, ( can be shifted to — £, 0, i.e., to 0, — £. 
Similarly, the hypergeometric equation has a group of F-homotopic automorphisms isomor- 
phic to (Z 2 ) 2 . These groups act on the 6 and 3-dimensional parameter spaces, respectively. 
For example, one of the latter actions is (a, b; c) 1— > (c — a, c — b; c), which is induced by an 
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F-homotopy at z = 1. From this F-homotopy follows Euler's transformation 2l|, §2.2] 

2 F x {a, b;c;z) = (l- z) c - a -\F 1 (c - a, c - b; c; z), (4.4) 



which holds because 2F1 is a local solution at z = 0, rather than at 2; = 1. If the singu- 
lar points of the differential equation are arbitrarily placed, transforming it to the Heun 
or hypergeometric equation will require a Mobius (i.e., projective linear or homographic) 
transformation, which repositions the singular points to the standard locations. A unique 
Mobius transformation maps any three distinct points in CP 1 to any other three; but the 
same is not true of four points, which is why (f)) has the singular point a as a free parameter. 



B. The cross-ratio orbit 



The characterization of Heun equations that can be reduced to the hypergeometric equa- 
tion will employ the cross-ratio orbit of {0, 1, d, 00}, defined as follows. If A, B,C, D e CP 1 
are distinct, their cross-ratio is 

iA B] °> D) = %-A)^-B) E ^ { °' °° } ' (45) 

which is invariant under Mobius transformations. Permuting A, B, C, D yields an action of 
the symmetric group 5 4 on CP 1 \ {0,l,oo}. The cross-ratio is invariant under interchange 
of A, B and C, D, and also under simultaneous interchange of the two points in each pair. 
Thus, each orbit contains no more than 4!/4 = 6 cross-ratios. The possible actions of S± 
on s E CP 1 \ {0, 1, 00} are generated by s 1— > 1 — s and s 1— > 1/s, and the orbit of s comprises 

8, 1-S, 1/8, 1/(1-8), S/(S-1), (S-1)/S, (4.6) 

which may not be distinct. This is called the cross-ratio orbit of s; or, if s = (A, B;C, D), 
the cross-ratio orbit of the unordered set {A, B, C, D} C CP 1 . Two sets of distinct points 
{Ai, Bi, Ci, Di} (i = 1, 2) have the same cross-ratio orbit iff they are related by a Mobius 
transformation. 
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C. Reminder of some of the properties of Heun's function 

Our aim will be to find an integral representation of the Heun function as a Frobenius' 
solution of the Heun equation, given in another form as follows js]: 

z(z -l)(z- a)y"{z) + {^{z - \){z - a) + 5z(z - a) + ez(z - l)}y'(z) 
+(a/3z-q)y(z) = 0, (4.7) 

The Frobenius' solution, noted Hl(a,q;a,f3,j,8;z) is the entire solution defined for the 
exponent zero at the point z — 0. It admits the power series expansion 



Hl(a, q; a, (3, 7, 5; z) = ^ c nZ n , 



n=0 



with \z\ < 1 and Cq = 1, C\ — and 7 7^ 0, —1, —2, 
The recursion relation is as follows: 



(4.8) 



a(n + 2)(n + 1 + l)c n+2 

= q + (n + l)(a + - 5 + (7 + 5 - l)a) + (n + l) 2 (a + 1) 
-in + a)(n + (3)c n = n > 0. 

The function Hl(a, q; a, /3, 7, 5; z) is normalised with the relation 

Hl(a, q; a, j3, 7, 5; 0) = 1. 



-n+l 



It admits the following important particular cases (|8J], p9, formula(1.3.9)): 

#7(1, a/3; a, /3, 7, 5; z) = 2 iq(a,/3,7; z) V<5 e C 
777(0, 0; a, /?, 7, 5; 2) = 2 7q(a, /3, a + (3 - 5 + 1; 2) V7 G C 
777 (a, a, /?, 7, a + /3 - 7 + l;z) = 2-^1(0;, /3, 7; -2), 



where 2-^1(0,/?, 7; 2) is the usual notation for the Gauss hypergeometric function. 



(4.9) 



(4.10) 



(4.11) 



D. Application of Heun's equation to our problem 



Finally we come to the punch line, why do we need these all and how does the Heun 
equation indeed find an application to our problem? The answer to this goes along the 
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following line: On us.ng m the difeentia. Rations for K and the equations (3.14) 
and (3.15) of |6j, acquire the form 

- {A + C)6 + k'+ -(A - C)k + mvr = 0, 

(C - A)U0 + 7T 7 + -(.A + C)7T + I71K = 0, (4.12) 

where and u are given in ( 13. 4ft . Clearly one can solve algebraically the second equation for 
ft. By differentiating the result one obtains also in terms of n, n' and 7r", and substitution 
of these in the first equation yields a second order ODE for 7r(/x). Now let us look at the 
system (14.121) in case of de Sitter spacetime. On inserting A and C from (15. ip below and 
passing to the globally defined variable z = cos 2 ^ (see Sec. Ill), we obtain the following 
differential equation for 7r: 



r, d r, d „ 

dz z dz 



n = 0, (4.13) 



where Pi in (14.131) is a quartic polynomial in z, i.e. 

P 2 = 4 [m 2 P 2 + l] z A - 4(2m 2 P 2 + 3)z 3 + 4(m 2 P 2 + 2)z 2 . (4.14) 

Similarly, P\ in (14.131) is a cubic polynomial in z, 

P a = 16 [m 2 P 2 + 1] 2 3 - 12 [2m 2 P 2 + 5] z 2 + 8 (m 2 P 2 + 2) z. (4.15) 

Last, P m (14.131) is a quadratic polynomial in z, i.e. 

P = (4m 4 - 19m 2 + 32m 2 P 2 + 9) 2 2 -(4m 4 - 14m 2 + 32m 2 P 2 + 2l) ^-3m 2 P 2 -6. (4.16) 

On making the substitution ir(z) = ^tt(z), (14.131) becomes an equation of the type 

z(z - l)(z - a)y"(z) + {(b + c+ l)^ 2 - [b + c + 1 + a(d + e) - e] z + ad}y'(z) 

+ {bcz-q)y{z) = 0. (4.17) 

Written in canonical form it reads as follows: 

d2y + (I + e + (b + c+1) - (d + e) \ dy + bcz - q ^ = Q 



dz 2 \z z— 1 2 — a / z(z — l)(z — a) 
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where the parameters in (14.181) take the values 

(m 2 R 2 + 2) 



a 



(m 2 R 2 + 1)' 
b = 2 + imR, 

c = 2 — imR, 

d = e = 3, 

(m 4 # 4 + 7m 2 R 2 + 10) 
q ~ (mm 2 + 1) ' { } 

The equation (14.171) is known as Heun's differential equation 0, Its solutions, here 
denoted by Hl(a, q; b, c, d, e; z), have in general four singular points as we said before, i.e. 
Zq = 0, 1, a, oo. Near each singularity the function behaves as a combination of two terms 
that are powers of (z — Zq) with the following exponents: {0, 1 — d} for Zq = 0, {0, 1 — e} for 
Zq — 1, {0, d + e — b — c} for z = a, and {b, c} (that is, z~ b or z~ c ) for z — » oo. 

We now insert into the second of Eq. (4.12) the first of Eq. (3.4), finding eventually 



K 



r'{ (A-C)((A + C) 2 -m 2 )-^(A + C)(m 2 + R- 2 ) n - (m 2 + R~ 2 )A , (4.20) 



where 

f = m{m 2 + R- 2 + 2C(C-A)), (4.21) 

and 7r and tt' are meant to be expressed through the Heun function Hl(a, q; b, c, d, e; z). 
Eventually, we will show in the next section that all weight functions can be therefore 
expressed through such Heun function. The material covered in the present section and in 
the previous two is not new, and most of it is appropriate only for a physics-oriented choice 
of four-dimensional de Sitter space. 



V. DICTIONARY OF WEIGHT FUNCTIONS FOR THE GRAVITINO PROPA- 
GATOR 

Here we will explicitly list all the weight functions as functions of z = cos 2 in order 
to analyze their qualitative behavior as a function of z and de Sitter radius R in the next 
section. Let us recall a few definitions in de Sitter space, where A and C are functions of 



the geodesic distance \i and are given by 11] 



A=— cot and C = — — ; — r--r, (5.1) 

R R Rsi-a^fR)' v ' 
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Since all other weight functions a, (3, 7, 5, e, 8, r, uo can be written in terms of the (tt, k) pair, 
and in the last section we have seen k can also be expressed in a form like f)4.20p . it is evident 
that all other 9 weight functions including k, i.e. a, /3, 7, 5, e, 9, r, u, k can be expressed in 
terms of 7r(/i) and 7r'(/x) only. 

We can also express 7r as a function of z and R only as 7r = ir(z) = 7r(/i = ±2i?cos~ 1 y / i). 
Similarly, by using a few of the familiar trigonometric identities, one can transform 7r'(/x) as 

n>(j J ) = ^>/z\T=z)ir>(z). (5.2) 

One can also write down the expressions of (A + C) and (A — C) in terms of z and R only 
as follows: 



*- c =Wi^- (5 - 3) 

Another function appearing quite frequently in our evaluation of all the weight functions is 
/, which can be also expressed as a function of z and R only as follows: 

f = m(m 2 + R- 2 + R~ 2 (l- z)- 1 ). (5.4) 

Now we start by listing all the weight functions in terms of tt(z) and tt'(z), bearing in mind 
that 

tt(z) — Hl(a, g; b, c, d, e; z), (5.5) 

7r(z) = \fz Hl(a, q; b, c, d, e; z), (5.6) 

where Hl(a, g; b, c, d, e; z) is the Heun function with arguments as defined before. One has 
therefore the lengthy formulae for all other weight functions written down in Eqs. (Al)-(A8) 
of the appendix. 



VI. QUALITATIVE BEHAVIORS OF THE WEIGHT FUNCTIONS 

Now using the series expansion (14.81) defined before one can numerically study the behav- 
ior of each weight function, by taking the first 10 terms of the infinite series (4.8). Indeed, 
dealing with an infinite number of terms is impossible, and one has therefore to resort to 
approximations, by truncating such a series. On taking less than 10 terms, we have found 
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minor departures from the pattern outlined below in figures 1 to 9, whereas on taking 15 
terms, the pattern in such figures is essentially confirmed. 

We draw for example all these weight functions in a two-dimensional plot vs z, in the 
range (0,1). The plots, which also include 7t(z), are as follows. 




FIG. 1: Two-dimensional plot of the weight function a{z). The curve has two branches, depending 
on whether one takes the + or — sign in (Al). One branch of a cuts the horizontal z-axis at 
the points z = 0.25,0.75, whereas the other branch of a cuts the horizontal axis at the points 
z = 0.12,0.6,0.92. Both branches approach the vertical axis, the first one cuts it near the value 
0.5, while the other has a vertical asymptote at z = 0.02. The two branches intersect each other 
at z = 0.15,0.69, 1; at these points the function a(z) becomes single- valued. 



i - 



-i - 




FIG. 2: Two-dimensional plot of the weight function (3{z). The curve has two branches, depending 
on whether one takes the + or — sign in (A2). One branch of (3 cuts the horizontal z-axis at 
the points z = 0,0.7,1, whereas the other branch of j3 cuts the horizontal axis at the points 
z = 0.18,0.67,1. The first branch never cuts the vertical axis, while the other has a vertical 
asymptote at z = 0.05. The two branches intersect each other at z = 0.15,0.69, 1; at these points 
the function (3(z) becomes single-valued. 
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FIG. 3: Two-dimensional plot of the weight function 7(2). The curve has two branches, depending 
on whether one takes the + or — sign in (A3). The first branch of 7 cuts the horizontal z-axis 
at the points z = 0.43,0.92 and the vertical axis at 1, and then it has a vertical asymptote at 
z = 0.95. The second branch of 7 cuts the horizontal axis at the points z = 0.47,0.76 and the 
vertical axis at 0.67, and then it has a vertical asymptote at z = 0.87. The two branches intersect 
each other at z = 0.15,0.5,0.69; at these points the function 7(2) becomes single- valued. The first 
branch of 7 has an absolute minimum, of negative sign, at z = 0.85. 



FIG. 4: Two-dimensional plot of the weight function 5(z). The first branch cuts the horizontal 
axis at z = 0,0.47,0.75, and the second branch cuts the horizontal axis at z = 0.45, 1. While the 
first branch never cuts the vertical axis, the second one cuts it at —2. The two branches intersect 
each other at z = 0.15,0.5,0.69, where 5 becomes single-valued. The first branch has a vertical 
asymptote at z = 0.85, whereas the second one does have the same at z = 1. 

As one can see, for values of z < 0.1, the main contribution to the gravitino propagator 
results from the weight functions a(z),/3(z), t(z) = u(z), whereas the other weight functions 
are sub-dominating. By contrast, when z e]0.8, 1[, the dominating contribution to the 
gravitino propagator results from the weight functions j(z),5(z),£(z),9(z),ir(z), while the 



-4 - 



-2 




15 



FIG. 5: Two-dimensional plot of the weight function s(z). The first branch cuts the horizontal 
axis at z = 0.2, 0.65, 0.93 and cuts the vertical axis at —1. The curve has an absolute minimum, of 
negative sign, at z = 0.85, and then reaches a vertical asymptote at z = 0.97. The second branch 
cuts the z-axis at z = 0,0.72. The two branches intersect each other at z = 0.15,0.69 and at these 
points e is a single-valued function. The first branch has a vertical asymptote at z = 0.95 whereas 
the second does the same for z in between 0.85 and 0.9. 



FIG. 6: Two-dimensional plot of the weight function 9(z). The first branch cuts the horizontal 
axis at z = 0,0.47,0.75 and the second cuts the horizontal axis at z = 0.42, 1. The first one never 
cuts the vertical axis, whereas the second one does it at the functional value 0.32. The second 
branch has a more pronounced absolute minimum, of negative sign, at z = 0.95. The two branches 
intersect each other at z = 0.15,0.5,0.69, where 6 is single- valued. The first branch has a vertical 
asymptote for z in between 0.85 and 0.9, whereas the second one does the same at z = 1. 

others remain sub-dominating. 
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FIG. 7: Two-dimensional plot of the weight function t(z) = u(z). The first branch cuts the 
horizontal z-axis at the points z = 0,0.7,1, and it never touches the vertical axis, whereas the 
second branch cuts the z-axis at z = 0.18, 0.67, 1, and it reaches a vertical asymptote at z = 0.05. 
The two branches intersect each other at z = 0.15,0.69, 1, where t(z) becomes single-valued. 



FIG. 8: Two-dimensional plot of the weight functions n(z) and tt{z). The n(z) curve passes through 
the origin and cuts the horizontal axis at z = 0.45, 0.83. The tt curve never passes through the 
origin, it cuts the vertical axis at the functional value 1 and it cuts the z-axis at z = 0.45, 0.83, 
where it also intersects the tt(z) curve. Beyond the point z = 0.8 the ir and tt curves become 
virtually indistinguishable. At z = 1 they both have a vertical asymptote. 

VII. CONCLUDING REMARKS 

Our paper has obtained the complete analytic structure of massive gravitino propagators 
in de Sitter space. In Sec. VI we have plotted all weight functions a, j3, 7, 5, e,9,r — u, n, k 
occurring in the gravitino propagator (jointly with fr) as a function of z in a two-dimensional 
plot where z = cos 2 (/i/2_R), ji being the geodesic distance between the points x and x', and 
R is the de Sitter radius. Although the series (4.8) has been truncated, it remains true 
that Sec. VI is the first attempt to display a supersymmetric propagator in de Sitter via 
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FIG. 9: Two-dimensional plot of the weight function k(z). The first branch cuts the z-axis at 
z = 0.29,0.7,0.97, while the second one cuts the z-axis at z = 0.68, 1. The first branch intersects 
the vertical axis at the functional value 1.35, whereas the second one does the same at the functional 
value 0.67. The two branches intersect each other at z = 0.15, 0.69, where n becomes single- valued. 

Heun functions. As we already said in Sec. I, further interest, from the point of view of 
mathematical methods, arises from the possibility to expand Heun functions in terms of 
hypergeometric functions As we said before, direct implications of our findings on the 
current understanding of the propagation of gravitinos in de Sitter space are as follows: there 
exist two ranges of values of z in which the weight functions can be divided into dominating 
and sub-dominating family. In other words, when z is smaller than 0.1, the weight functions 
a,/3,T = u are dominating while the others are sub-dominating. By contrast, when z is 
very close to 1, the weight functions 7, 5, e, 9, tt take much larger values. 

The plot range is between and 1 for z, which is indeed the only admissible region, since 
the squared cos function lies always between and 1. Note that the plot of tt is basically 
nothing but the plot of the Heun function with properly defined coefficients, and the plot of 
7r is y/z times the Heun function. 

The numerical analysis of Sec. VI, as we already said therein, has been performed by 
taking only the first 10 terms of the infinite series representing the Heun function, by ap- 
plying the Frobenius' method. If one goes on by taking more terms, one can get even more 
accurate results, but roughly the qualitative features remain the same. The task of plotting 
Heun functions is technical but not easy, since the modern computer packages still run into 
difficulties. Thus, our efforts can be viewed as preparing the ground for a more systematic 
use of Heun functions in fundamental theoretical physics. The flat-space limit is instead a 
considerable simplification, since the functions A and C in (5.1) are then found to reduce to 
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A = C = — ^, and the formulae in the appendix are therefore considerably simplified. 

It also remains to be seen whether the familiarity acquired with Heun functions will 
prove useful in studying gravitino propagators in other backgrounds relevant for modern 
high energy physics. 

APPENDIX A: EXPLICIT FORM OF THE WEIGHT FUNCTIONS 

The weight functions obtained in Sec. V read, explicitly, 



a(z) = -2mC(m 2 + R~ )f (z) x 



1 



(A - C)((A + C) 2 - m 2 ) - -(A + C)(m 2 + R~ 2 ) 
-((m 2 + R- 2 )[(A + C) 2 - m 2 ] - 4)tt(z), 



tt(z) ± (m 2 + R' 2 ) ^ Z{ \ Z K '(z) 



R 



(Al) 



(5{z) = AmC(m 2 + iT 2 )" 1 /" 1 ^) x 

(A - C)({A + C) 2 - m 2 ) - \{A + C)(m 2 + R~ 2 ) 
+2(m 2 + R- 2 )- 1 [(A + C) 2 - m 2 ] n(z), 



n(z)±(m 2 + R- 2 )^- Z) -> 

R 



A*) 



(A2) 



7 (*) = -{m 2 + R- 2 )'' [(A - C) 2 - m 2 } f~\z) x 
(A - C)((A + C) 2 - m 2 ) - \{A + C)(m 2 + R7 2 ) 

Z 

-2mC{m 2 + R~ 2 y\(z), 



ir(z) ± (m 2 + R~ 



R 



-7T (Z) 



(A3) 



8{z) = -6(m 2 + R- 2 )- 1 [(A - C) 2 - m 2 ] f~h 



z) x 



1 



(A - C)((A + C) 2 - m 2 ) - -(A + C)(m 2 + R~ 2 ) 

z 



+ l2mC(m 2 + R- 2 )- L ir(z), 



n(z) ± (m 2 + R- 2 ) 



R 



TX (Z 



(A4) 



e( z ) = -(m 2 + R- 2 ) 



2\-l 



(A-C) 2 + - + m 2 



r\z) x 



{A - C){{A + C) 2 - m 2 ) - -(A + C)(m 2 + R~ 2 ) 

z 



+2mC(m 2 + R- 2 y L ir(z), 



Tr(z)±(m 2 + R- 2 )y /z{1 ~~ ] 



R 



n'(z) 



(A5) 
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B{z) = (m 2 + R^Y 1 [{A - Cf - m 2 } f~\z) x 

(A - C){{A + Cf - m 2 ) - Ua + C)(m 2 + iT 2 ) 
-2mC{m 2 + R~ 2 y 1 7r{z), 



n(z)±(m 2 + R- 2 )^- Z) ~> 

R 



7T (Z 



(A6) 



t(z) = 2mC(m 2 + R~ z )- f ( z ) x 

(A - C)((A + C) 2 - m 2 ) - ^{A + C)(m 2 + R~ 2 ) 
+ (m 2 + R- 2 Y l [(A + C) 2 - m 2 ] ir(z) = lu(z) } 



ir(z)±(m +R )- — 

R 



A*) 



(A7) 



K (z) = r\ Z ) x 



1 



(A - C)((A + C) 2 - m 2 ) - -{A + C)(m 2 + R~ 2 ) 

_ 2 y/z(T=z), 



ir(z) 



(A8) 



These exhaust all the weight functions multiplying the invariant structure present in the 
gravitino propagator, written explicitly in terms of a Heun function and its derivative. 
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